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Abstract 

The formation of a polaron quasiparticle from a bare electron is studied in the framework of 
the Holstein model of electron-phonon coupling. Using Schrodinger's formalism, we calculate the 
time evolution of the distribution of the electron and phonon density, lattice deformation, and the 
electron-phonon (el-ph) correlation functions in real space. The quantum dynamical nature of the 
phonons is preserved. The polaron formation time is related to the dephasing time of the continuum 
of unbound phonon excited states in the spectral function, which depends on the phonon frequency 
and the el-ph coupling strength. As the el-ph coupling increases, qualitative changes in polaron 
formation occur when the one-phonon polaron bound state forms. In the adiabatic regime, we find 
that a potential barrier between the quasi- free and heavy polaron states exists in both 2D and 3D, 
which is crucial for polaron formation dynamics. We compare to recent experiments. 

PACS numbers: 71.38.-k, 63.20-e, and 78.47.+p 
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The polaron formation problem not only bears fundamental interest in understanding 
polarization in condensed matter, but also determines the electronic and optical properties 
in some applied materials, such as manganites |1] and conductive organic oligomers Q|. 
Recent advances in ultra-fast time-resolved spectroscopy have made it possible to investigate 
physical phenomena on the time scale of a molecular vibration or optical phonon in crystals. 
Important aspects of chemical reaction dynamics have been revealed by observation of the 
motion of the oscillating atoms [3|. Observations of the femtosecond dynamics of polaron 
formation have recently been reported jj, 0, (], B, Is]. In contrast, theoretical development 
in this subject is somewhat behind. Although the theoretical research on polaron physics 
spans six decades |9|, ll(|, the dynamics of lattice relaxation leading to a heavy polaron state 
remains far from being und e rstood . An adlabati c polaI on ^ Q suited that te e 
is a potential barrier between the delocalized and the self-trapped polaron states in 3D but 
not in ID or 2D. It is generally believed that the potential barrier tends to impede the 
polaron formation. However, more careful thought is needed to go beyond adiabatic limit 
(i.e. ujq > 0) where all eigenstates are delocalized. It is not clear whether the barrier slows 
down or forbids the formation of a polaron. 

Measuring polaron formation time has been the subject of the above recent experiments 



]. The time to form a polaron is found to be less than a picosecond in strongly 
coupled systems, on the order of a phonon period. Our goal is to provide a theory for the 
following fundamental questions: (1) How are phonon excitations triggered and how do they 
evolve into the correlated phonon cloud of the polaron quasiparticle? (2) How much time 
does it take to form a polaron? (3) What is the effect of dimension on the dynamics of 
polaron formation? 

To understand the formation process of a polaron, we examine how the bare particle wave 
function time evolves into a polaron quasiparticle. One approach is to construct a variational 
many-body Hilbert space including multiple phonon excitations, and to numerically integrate 
the many-body Schrodinger Equation, i^- = Hip in this space (l2|. Thus the full many- 
body wave function can be obtained at early times. The main approximation is the size of 
the variational space, which can be increased systematically until convergence is achieved. 
This method includes the full quantum dynamics. Alternative treatments, such as the 
semiclassical approximation Q], can be inaccurate when applied to the present problem. 
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We solve for the dynamics of the Holstein Hamiltonian 



H — H e i + H e i_ p h + H p h 
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where cj creates an electron and at creates a phonon on site j. The parameters are the 
nearest-neighbor hopping integral t, the el-ph coupling strength A, and the optical phonon 
frequency u . 

Figure E shows snapshots of polaron formation at weak coupling. An initial bare electron 



wave packet is launched to t 



re right as shown in panel (a). (This initial condition is relevant 



to the recent experiments |4 



and to electron injection from a time-resolved STM 



tip 14j]). In panel (b) the electron is not yet dressed and thus is moving roughly as fast as 
the free electron (dashed line). In addition, there exists a back-scattering current (which 
evolves into a left-moving polaron later) on the left side of the wave packet (green and 
black curves). In panel (c) after an elapsed time of one phonon period, the electron density 
consists of two peaks. The peak on the right (black arrow) is essentially a bare electron. 
The peak on the left is a polaron wave packet moving more slowly. The velocity operator 

is defined as Vj = , 2J jj+ 1 where j is the site index and J is the current operator 

Jjj+i = —i et (c]cj + i — H.c.) . (Vj) is shown in green. As time goes on, the bare electron 
peak decays and the polaron peak grows. Some phonons are left behind (blue line), mainly 
near the injection point. These phonons are of known phase with displacement shown 
in red. Some phonon excitations travel with the polaron (magenta). Finally, a coherent 
polaron wave packet is observed when the polaron separates from the uncorrelated phonon 
excitations. 

There are regimes where the polaron formation time is a calculable constant of order 
unity times a phonon period T , as seen in some experiments and shown in Fig. ^ but 
there are other regimes where the phonon period is not the relevant timescale. The limit 
hopping t — ^ is instructive. After a time To/4, the expectation of the phonon coordinate 
(Xj) on the electron site has the same value as a static polaron. It is tempting (but we 
would argue incorrect) to identify this as the polaron formation time. At later times, (Xj) 
overshoots by a factor of two, and after time To, (Xj) and all other correlations are what 
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they were at time zero when the bare electron was injected. The system oscillates forever. 
In general an electron emits phonons enroute to becoming a polaron, and we propose that 
the polaron formation time be defined as the time required for the polaron to physically 
separate from the radiated phonons. The polaron formation time for hopping t — > is thus 
infinite, because the electron is forever stuck on the same site as the radiated phonons. 

An electron injected at several times the phonon energy uq above the bottom of the 
band is another instructive example. The electron radiates successive phonons to reduce its 
kinetic energy to near the bottom of the band, and then forms a polaron. For weak el-ph 
coupling A, the rate for radiating the first phonon can be computed by Fermi's golden rule, 

! ^ 

T FGR = *- . • 11 \i (2) 

where kf is the electron momentum after emitting a phonon. The phonon emission time can 
be arbitrarily longer than the phonon period T for small A. For strong coupling, the rate 
approaches t^q = X/% because the spectral function smoothly spans numerous narrow bands 
and its standard deviation is equal to A. Our numerical results agree with the perturbation 
theory (not shown). 

We now consider polaron formation in more detail. After injecting a bare electron at 
time zero, the wavefunction at later times r is 

oo 

l^(r)) = Ee-^ T |^)(^|cl|0), (3) 

3=1 

where \^fj) are a complete set of total momentum k eigenstates of the system of one electron 
coupled to phonons. There are several distinct types of states contributing to the infinite 
sum: (A.l) The state \k) of a momentum k polaron, corresponding to the quasiparticle 
pole. (A. 2) The states \k — q; q) corresponding to a polaron of momentum k — q and an 
unbound phonon of momentum q, for any q. Similarly for two unbound phonons, etc. If 
the electron-phonon coupling is sufficiently strong, there are additional states in the sum: 
(B.l) A polaron excited state consisting of a polaron and an additional bound phonon of 
total momentum k 15j ]. designated This is a second type of (excited) quasiparticle 

pole that is also split off from the continuum. (B.2) The states \(k — q)^; q) corresponding 
to an excited state polaron of momentum k — q and an unbound phonon of momentum 
q for any q. Similarly for two unbound phonons, etc. For stronger el-ph coupling, more 
highly excited polaron states corresponding to bound states of a polaron and two or more 



additional phonons, \k^ n ^), . . . enter the sum Q]. The branching ratios into the various 
channels are calculated in Ref. jl?! ]. 

From Eq. |3J the amplitude to remain in the initial state after time r, (^(t) |c|.|0) , is given 
by the Fourier transform of the spectral function 

oo 

A{k^) = Y.\(*i\Am 2 K"-^)- (4) 

i=i 

The numerically determined spectral function at fairly weak coupling is shown in Fig. |21 
We compute the spectral function by a Lanczos algorithm. The sum rule N(uj) shows that 
about 10% of the total spectral weight is at energies beyond the range plotted. There is a 
quasiparticle pole corresponding to (A.l) above, and a group of states that is approaching 
an approximately Lorenzian continuum as the number of sites increases, corresponding to 
(A. 2) above. The coupling X/ujq = 0.8 is too weak to form bound quasiparticle excited 
states \k^'). If the spectral function were a pure Lorenzian, a measurement would yield a 
exponential decay of the initial state, with the polaron formation time t p the inverse width 
of the Lorenzian. Since the spectral function has a quasiparticle pole as well, an experiment 
would measure the probability to remain in the initial state 

P(t) = a \ + a\e- 2hT + 2a x a 2 e- hT cos[(cui - ui 2 )t\. (5) 

This form already shows some complications, with an additive constant, a pure exponential 
decay, and an exponential decay half as fast multiplied by a cosine oscillating at the energy 
difference between the quasiparticle pole and the center of the Lorenzian. Decaying oscilla- 
tions in polaron formation (actually the formally equivalent problem of an exciton coupled 
to phonons [18]) have been observed in a pump-probe experiment 0] which measures reflec- 
tivity after a bare exciton is created. The observed oscillatory reflectivity was interpreted 
as the lattice motion in the phonon-dressed (or self-trapped) exciton level. Assuming the 
modulation in the exciton level goes as AE = —Xc^CjXj where Xj is lattice displacement, 
the model Hamiltonian applies directly to the experiment. We calculate the corresponding 
el-ph correlation function x = ( c ] c j( a j + a ]))j which is shown in Fig. El The result is in 
satisfactory agreement with experiment p. (Numerical calculations in Fig. El thro ugh Fig.03 
are performed on an extended system, not a finite cluster.) 

Figure 0] shows the spectral function at stronger coupling than Fig. El A polaron and 
two excited state polaron poles, along with three continua are shown. There is additional 



structure at higher energy (not shown). The probability decay -P(t) for this spectrum is 
considerably more complicated, and includes oscillating terms that do not decay to zero (at 
zero temperature) from the ground and excited polaron poles beating against each other. 

Next we discuss the role of dimensionality. The effect of dimensionality on static proper- 
ties has been studied previously The eigenvalues of the low-lying states (with zero total 
momentum) are shown as functions of A in Fig. EJ It is found that the energy spectra in 
D > 1 are qualitatively different than in ID. The ID polaron ground state becomes heavy 
gradually as A increases. However, in 2D or 3D, the ground state crosses over to a heavy 
polaron state by a narrow avoided level crossing, which is consistent with the existence of 
a potential barrier In the inset of Fig. |SJ ipi and ^4 are nearly free electron states; 

ip2 and ip 3 are heavy polaron states. The inner product KV'ilV^)! is equal to 0.99. (A simi- 
lar avoided level crossing also occurs in 2D.) Just right of the crossing region the effective 
mass (approximately equal to the inverse of spectral weight) of the first excited state can 
be smaller than the ground state by 2 or 3 orders of magnitude, while their energies can 
differ by much less than ujq. As a result, there is no optical phonon of the correct energy for 
the light electron to emit and become a heavy polaron, and the polaron formation time at 
T = would be infinite unless low energy acoustic modes are included in the model. 

In summary, we have calculated the time evolution of the many-body wave function 
and find that the bare electron evolves into a polaron quasiparticle by emitting phonons. 
The excess kinetic energy is used to excite uncorrelated phonons. The question "how long 
does it take a polaron to form?" may not have a simple answer, given the potentially 
complicated form of P{t). This function has been calculated numerically, and at zero 
temperature depends on the parameters of the Hamiltonian, the spatial dimension, the 
initial bare electron momentum k, the final polaron momentum, and the possible existence 
of bound polaron excited states. A further complication is that decay out of the initial state 
need not be synonymous with decay into a polaron final state, as seen from the weak coupling 
A regime, Eq. El In addition, we show that a tunneling barrier between the quasi-free and 
heavy polaron state exists in both 2D and 3D (when < ujq <C 2Dt). As a consequence, the 
tunneling barrier inhibits the formation of a polaron in the crossover regime. Our approach 
is intuitive and can be extended to other types of quasiparticle formation, such as the 
vibrational exciton |_2Q;] and spin polaron. 
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FIG. 1: Snapshots of the polaron-formation process, for hopping t = Uq = 1, and A = 0.4. 
The calculation is performed on a 30-site periodic lattice. Time is measured in phonon periods. 
Black: electron density (cjcj); Blue: phonon density (ajctj); Red: lattice displacement (a,j + 
Green: velocity in units of lattice constant per phonon period; Magenta: el-ph correlation function 
(CjCjOjdj); Dashed: free-electron wave packet for reference. For clarity, the origins of the red 
and green curves are offset by 0.1 and their values are rescaled by a factor of 0.2 and 0.05/(2-71") 
respectively. 
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FIG. 2: Spectral function at weak coupling. r p is the polaron formation time. 
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FIG. 3: The on-site correlation function \ as a function of time. The initial energy of the bare 
electron (or exciton) is Ei = —0.7. When A = 1.0, the result most resembles the experimental data 
5]. For all curves, ojq = 0.5 and t = 1. 
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FIG. 4: Panel (a): spectral function at strong coupling. There exist three quasiparticle poles. 
Shaded areas (1) and (2) correspond to continuum states (A. 2) and (B.2) respectively, (b): Quan- 
tum beat formed by multiple quasiparticle poles. 
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FIG. 5: Eigenvalues of low-lying states as functions of coupling constant in ID through 3D. 
Hopping t = 1 in all panels. In the adiabatic regime in higher dimensions, the ground state (thick 
solid lines) shows a fairly abrupt change in slope. In the 3D panel, ipi and ip^ are a lightly-dressed 
quasiparticle state; ^2 and ^3 are a heavy polaron state. The dashed lines are the beginning of 
the lowest continua. 
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